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We calculate the tunneling conductance spectra of a normal metal / insulator / triplet superconductor from 
the reflection amplitudes using the Blonder- Tinkham-Klapwijk (BTK) formula. For the triplet superconductor 
we assume one special p-wave order parameter having line nodes and two two dimensional /-wave order param- 
eters with line nodes breaking the time-reversal symmetry. Also we examine nodeless pairing potentials. The 
tunneling peaks are due to the formation of bound states for each surface orientation at discrete quasiparticles 
trajectory angles. The tunneling spectra can be used to distinguish the possible candidate pairing states of the 
superconductor Sr2Ru04. 



1. Introduction 

The resent discovery of superconductivity in 
Sr2Ru04 has attracted much theoretical and ex- 
perimental interest Muon spin rotation ex- 
periments show that the time-reversal symmetry 
is broken for the superconductor S^RuO^ [||. 
Knight-shift measurements show no change when 
passing through the superconducting state and is 
a clear evidence for spin triplet pairing state || . 
Also specific heat measurements support the sce- 
nario of line nodes within the gap as in the high 
T c cuprate superconductors [Q. 

In tunneling experiments involving singlet su- 
perconductors both line nodes and time-reversal 
symmetry breaking can be detected from the V- 
like shape of the spectra and the splitting of the 
zero energy conductance peak (ZEP) at low tem- 
peratures respectively 

Also the properties of ferromagnet - insulator 
- superconductor with triplet pairing symmetry 
have been analysed j|. It is found that the 
bound states are suppressed and hence the tun- 
neling conductance peaks are eliminated with the 
increase of the exchange field. 

In this paper we will use the Bogoliubov-de 
Gennes (BdG) equations to calculate the tunnel- 
ing conductance of normal metal / triplet super- 
conductor contacts, with a barrier of arbitrary 
strength between them, in terms of the probabil- 
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ity amplitudes of Andreev and normal reflection. 
For the triplet superconductor we shall assume 
three possible pairing states of two dimensional 
order parameter, having line nodes within the 
Ru02 plane, which break the time-reversal sym- 
metry. The first two are the 2D /-wave states 
proposed by Hasegawa et al, having B\ g x E u 
and £?2g x E u symmetry respectively. The other 
one is called nodal p-wave state and has been pro- 
posed by Dahm et al [|lTj| . This pairing symmetry 
has nodes as in the Bi g x E u case. Also we will 
consider the nodeless p-wave pairing state pro- 
posed by K. Miyake and O. Narikiyo fl2]| . 



2. Theory for the tunneling effect 

We consider the normal metal / insulator / su- 
perconductor junction shown in Fig. |l|. The ge- 
ometry of the problem has the following limita- 
tions. The particles move in the xy-plane and the 
boundary between the normal metal (x < 0) and 
superconductor (x > 0) is the yz-plane at x = 0. 
The insulator is modeled by a delta function, lo- 
cated at x = 0, of the form VS(x). The tempera- 
ture is fixed to K. We take the pair potential as a 
step function i.e. A S s(k, r) = <d(x)A S s(9), where 
k x ,k y = cos 9, sin 9, and s,s are spin indices. 

Suppose that an electron is incident from the 
normal metal to the insulator with an angle 9. 
The electron like (hole) like quasiparticle will ex- 
perience different pair potentials A S s(9)(A S s(tt — 
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Figure 1 . The geometry of the normal metal / 
insulator / superconductor interface. The vertical 
line along the y-axis represents the insulator. The 
arrows illustrate the transmission and reflection 
processes at the interface. 9 is the angle of the 
incident electron beam and the normal. 



9)). The amplitudes for Andreev and normal re- 
flection are obtained by solving the BdG equa- 
tions. Using the matching conditions of the wave 
function at x = 0, */(0) = *//(0) and *^(0) - 
^(0) = (2mV/h 2 )^ I (0), the Andreev and nor- 



mal reflection amplitudes R a 
obtained 



\a\ 2 ,R b = \bf 



Ra 



Rb — 



a%\n+\ 2 



|1 + (o-jv - l)n + n-4>-<p* + \ 2 ' 
(f - <t n )\1 - n+n^(j)^(j)* + \ 2 



The BCS coherence factors are given by 



= [l + VE 2 -\A ± (0)\yE]/2, 



[1 - VE 2 - |A±(0)| 2 /£]/2, 



(1) 
(2) 

(3) 
(4) 



and n± — v±/u±. The internal phase com- 
ing from the energy gap is given by <f>± = 
[A±(0)/|A±(0)|], where A+(0) = A(9) (A_(9) = 
A(7r — 9)), is the pair potential experienced by 
the transmitted electron-like (hole-like) quasipar- 
ticle. A(0) = A Ti (6>) = A iT (0), since the Cooper 
pairs have zero spin projection i.e. d || z. 



The tunneling conductance, normalized by that 
in the normal state is given by H 



a(E) 



j:i 2 /2 d8(ai(E,8)+<T l (E,6)) 
Jli 2 /2 d92a N 



(5) 



According to the BTK formula the conductance 
of the junction a s (E,9), for spin s =t,|, is ex- 
pressed in terms of the probability amplitudes a, 
and b as 



W S (E,9) = 1 + R a -R b . 



(6) 



The transparency of the junction er/v is connected 
to the barrier height V by the relation 



4 cos 2 



<7n = 



Z 2 + 4 cos 2 9 ' 



(7) 



where Z — 2mV/h Uf, denotes the strength of 
the barrier. 

The pairing potential is described by a 2 x 2 
form 



A s ^(k) 



-d x (k) + id y (k) 4(k) 

4(k) d x (k) +id y (k) 



in terms of the d(k) = (d x (k), dj/(k), <i z (k)) vec- 
tor. 

We consider the following pairing symmetries 
for Sr 2 Ru0 4 . 

a) In the first 2D /-wave state B\ g x E u d z (k) = 
^oik 2 . — k 2 )(k x + iky). This state has nodes at 
the same points as in the d x 2_ y 2-wave case. 

b) For the second 2D /-wave state Bi g x E u 
d z Qc) = Aok x k y (k x + iky). This state has nodes 
at 0, 7t/2, 7t, 37t/2 and has also been studied by 
Graf and Balatsky §J|. 

c) In case of a nodal p-wave superconductor 
d z Qs) = j- 1 [sm.(k x a,) + «sin(fc y a)], with k x a = 
i?7rcos(0 — (3) and k y a = i?7rsin(0 — (3), sm = 
\/2sin -3= = 1.125, and R = 1 in order to have a 
node in A(0) jy|. This state has nodes as in the 
B2g x E u state. 

d) In case of a nodeless p-wave superconduc- 
tor, proposed by K. Miyake, and O. Narikiyo Jl2] ] 
d z (k) = Aq. [sin(fca;a) + ism(k y a)], with k x a = 
Rircos(9 — (3) and k y a = i?7rsin(0 — f3), sm = 
V2sin^ = 1.125, and R = 0.9. This state does 
not have nodes. 



(8) 
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Figure 2. Normalized tunneling conductance 
a(E) as a function of E/Aq for different orien- 
tations (3 = (solid line), 7r/8 (dotted line), 7r/4 
(dashed line). The temperature is T = 0, Z = 10. 
The pairing symmetry of the superconductor is 
(a) Big x E u , (b) B 2g x J5 U , (c) nodal p-wave, (d) 
nodeless p-wave. 



3. Results 

In Fig. |^ we plot the tunneling conductance 
cr(E) as a function of E/A for Z = 10, for dif- 
ferent orientations (a) (3 = 0, (b) ir/8, (c) 7r/4. 
The pairing symmetry of the superconductor is 
-Bi s x E u in Fig. |2|a, i?2 9 x E u in Fig. ||b, nodal 
p-wave, in Fig. gp, nodeless p-wave in Fig. ^d. 
The temperature is fixed to OK. In all the cases 
we see the presence of a large residual density of 
states within the energy gap and peaks due to the 
sign change of the pair potential at discreet val- 
ues of 9, for fixed /3. Also the linear increase with 
E is consistent with the presence on line nodes in 
the pairing potential. 

Generally the spectra for the three pairing sym- 
metries with line nodes depends strongly on the 
position of the nodes in the pairing potential and 
the orientation angle (3. The spectra for angle 
(3 in the Bi g x E u seen in Fig. ||a case is iden- 
tical to the spectra of B 2g x E u in Fig. for 
angle 7r/4 — j3, since the node positions for the 



Figure 3. Bound state energy E p (in units of 
Aq) for (3 = 0,7r/8, 7r/4 as a function of 9. The 
temperature is T = 0. The pairing symmetry of 
the superconductor is (a) B ig x E u , (b) B 2g x E u , 
(c) nodal p-wave, (d) nodeless p-wave. 



two symmetries differ by 7r/4. The nodal p-wave 
case in Fig. ||c has the same nodal structure as 
the B 2g x E u case and we see that the spectra 
for these two candidates are similar. For nodeless 
pairing states a subdap or a full gap opens in the 
tunneling spectra. This is seen in Fig. |^d, for the 
nodeless p-wave pairing state. 

The peaks in the tunneling conductance are ex- 
plained from the formation of bound states within 
the gap. The bound state energies Ep , are given 
from the values of E where the denominator of 
Eqs. (|l],||) vanishes. In this case the Andreev re- 
flection coefficient is equal to unit, and the effect 
of the boundary is turned off. The corresponding 
equation is written as |tJ 

^* + n+n_\ E=Ep = 1.0. (9) 

Bound states are formed because the transmitted 
electron-like and hole-like quasiparticles feel dif- 
ferent sign of the pairing potential. These bound 
states occur for a given orientation /? at discreet 
values of 9, as seen in Fig. [| where the bound 
state energy E p is plotted for f3 = 0, 7r/8,7r/4 as 
a function of 9 for the various pairing states. For 
a given value of (3, the conductance <r(E, 9) = 2 
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at the angles 9 where bound state occurs and the 
tunneling conductance cr(E) is enhanced due to 
the normal state conductance cr/v- The residual 
values of the tunneling conductance within the 
gap is due to the bound states and the peaks are 
formed at energies where the number of bound 
states is increased. As seen in Fig. | (a) for the 
pairing state B\ g x E u for j3 = 0, at the energy in 
the interval 0.27 < E < 1 only one bound state 
occurs. At E — 0.27 two more bound states are 
formed and the peak seen in Fig. || (a) is due 
to these new bound states. The same happens 
for P — 7r/8 at E — 0.1 and E = 0.5, where an 
increased number of bound states is formed, as 
seen m Fig | (a) (dotted line). Also for (3 = tt/4 
in Fig U (a) (dashed line) the conductance peak 
is formed for energy close to E = 0.8, where two 
bound states are formed. In the Big x E u pair- 
ing state the bound states and also the position 
of the peaks for an angle (3 is the same as in the 
Bi g x E u pairing state for the angle 7r/4— f3. In the 
nodal p-wave state seen in Fig. || (c) the bound 
states are symmetric to the B 2g x E u case with re- 
spect to 8 = 0. However this does not influences 
the energy levels which occur almost at the same 
position as in the Bi g x E u -w&ve case. In the 
nodeless p-wave pairing state, for j3 = w/8 close 
to E = we have a pair of new bound states and 
the peak seen in Fig. [2] (d) is due to these new 
bound states. 

We calculated the tunneling conductance in 
normal metal / insulator / triplet superconductor 
using the BTK formalism. We assumed nodal and 
nodeless pairing potentials, breaking the time- 
reversal symmetry. The residual values are due 
to the formation of bound states for each /3 at 
discreet values of the angle 9 for energies within 
the gap. The peaks occurs at the energies where 
an increasing number of bound states is formed. 
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